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In our subgroup, we focused on penguin physics from various different angles. Our 
discussions included: (1) A method to extract one of the angles of the unitarity 
triangle </>2; (2) Methods to extract <j):j from B — » Kir decays; (3) Effects of non- 
minimal SUSY on B and K decays; (4) Understanding large branching ratios for 
B — > + K(K*) decays; (5) New calculation for hadronic matrix elements which 
are needed to compute — . 



1 Introduction 

Our session was given a rather general title: "Weak Decays, CKM and CP 
Violation" . It is a big field and we can not do justice to any of these subjects 
if we try to cover everything. For this reason, we decided to concentrate our 
discussion on penguin physics. 

Year 1998 was a very good year for penguin physics - year 1999 promises 
to be even better for flavor physics in general. 

1. We ara^upplied with branching ratios on hadronic two body decays from 
CLEold: 
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2. New results on ^ was promised, and in fact, just after the meeting the 
result from E83213 was announced. The result is considerably larger than 
that of previous Fermilab result!: 



Re- = 



(23 ± 6.5) X 10"'' NA31 
e ~ \ (28.0 ± O.SOstat. ± 0.26,j,,t. ± O-lOucstat.) x 10-^ E832 

This result is now in agreement with that from CERN@ and establishes 
non-vanishing direct CP violation. 

3. Belle and Babar collaborations should start taking data on much antici- 
pated large CP violation in B decays. Along the way, there will get lots 
of data on B decays. 

4. A new K meson factory at Da$ne should start taking data this year. 

Who knows, by year 2001, we may have a positive signal for New Physics. Much 
of above experimental development demands better theoretical understanding 
of penguins. 

2 How big are penguins? 

Let us illustrate the importance of penguins by giving a hand- waving argument 
based on the experimental result Eq. ([^). Less hand- waving argument is 
presented in Gronau's contributionQ. 

The amplitudes for tree and penguin contributions for Ki: decay mode 

are: 

P{Knl = ^V:,V,4Cl{^i)Q';,iK^T) + C2(A*)g^2(i^^)] 
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P{Kn)t = -^{-V;,Vts)Y.ai^^)Qs^{K7^). (2) 

^ i=3 

To simplify our notation, set: 

TiKn) = V:,V„,r; P{K7r), = V*,V,,P^; P{KTT)t = -V*,Vt,P*. (3) 



If penguin diagram gave negligible contribution, the entire two body decays 
occur through Fig. m(a). Then B TTTT decay goes through a diagram in which 
the s quark in Fig.l^(a) is replaced by a c? quark. For a rough argument, we 
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Figure 1: Quark diagrams for Jf+Tr'^ and tt+tt" decays, (a) The tree graph contribu- 

tion, (b) Penguin contribution. 



ignore SU(3) breaking in the hadronic matrix elements. Then, B tttt is 
given by T(t:ti) ^ X^T. Since T{Kii) IT{-k'k) ^ A, we would expect: 

Br(B ^ i\:7r)/Br(B ^ tttt) - 0(A2). (4) 

Experimentally this is not so. This indicates that the P{KTr) amplitude is at 
least as large as the T{-ktt). 

For B Kit, P{Ktt) is ^(A^) and T{Kn) is ©(A"*). So, P{KTr) must be 
a major contributor to the decay amplitude. If P{Ktt) ~ T(7r7r), then 

T>c I T>t 

--OiX), (5) 



r 

the penguin contribution is considerably larger than what a naive estimate of 
the loop graph would suggest: 

-— T^T^log O 0.01 . 6 
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Figure 2; Angles of the unitarity triangle are related to the phases of the KM matrix. The 
right hand rule gives the positive direction of the angle between two vectors. 



Gronau went through less hand-waving analysis and obtained 



0.3 ±0.1. (7) 



Large penguin contribution is not always welcome. For example, they play 
a role in so called "penguin pollution" which causes hadronic uncertainty in 
determining 02 and 03 . For notations see Fig. Problem penguins cause, 
however, does not compare with richness they bring to flavor physics. Unlike 
in K decays where effects of tree graphs dominate, in B physics, quantum loop 
effects via penguins is often a leading contribution. This gives us a window of 
opportunity to look for effects beyond the standard model - as they arc likely 
to contribute through loop effects. Anticipating this possibility, we had the 
following discussions: 

1. Reviews of penguins in B decays by M. Gronau. 

2. New remarks on the determination of (pi and 02 by L. Oliver. 

3. A critical look at 03 from B Kit by R. Fleischer. 

4. Model independent anlysis of i? ^ Ktt decays and bounds on the weak 
phase 03 by M. Neubert. 

5. Analysis of B -> r]K{K*) and B rj'K{K*) by D-D. Du. 

6. Effects of SUSY particles in B and K decays by A. Masiero. 

7. Chiral methods and predictions for K ^ tttt hy E. Paschos. 



Here we use the notation which was introduced when the unitarity triangle was first 
discussed in the context we use today B. 
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3 Taming the penguins 



How to get around the penguin pollution and extract the value of <j)2 has been 
reviewed by Gronau. Oliver has presented an alternative approach which may 
be useful. In his approach, (j)2 is expressed as a function of experimentally 
measurable quantities in _B — s- tttt decay, plus one other parameter. It can be, 
for example ^, obtained from B Kn, mentioned above. 

The time dependent B(B) tt+tt^ decay rates are given by: 



r ( ^ ^ TT+TT- ) (t) 



\A\' + \AfT{\Af-\Af) cos{AMbJ) 

q- 



± 2Im ( - AA* ] 8in{AMBj) 



(8) 



There are three experimental observables: 



\Ai Ml, and Im ^AA' 



Theoretically, we can write 



1 



(9) 



Here M" = T" - 7"=, and M* ^ - V. These are related to previously 
introduced matrix elements except for the SU(3) breaking corrections. 
Theoretical unknowns are 



|M"| 



arg 



Since there are 4 theoretical parameters and only 3 experimental observables, 



> Kit 



we can not solve for 02- We can, however solve for 02 as a function of, e.g. 

We can, most likely obtain this parameter from Eq. ( |l6| ) looking at B 
decays. Further study is necessary to see how the error from SU(3) symmetry 
breaking will affect the determination of (f>2- Also, there are some ambiguities 
coming from the sign oLa square root as well as that coming from 02 ± tt. For 
details see Oliver's talkH. 
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4 Getting maiximum out of B ^ Kir, tttt decays 

The CLEO coUaboration has recently reported the observation of _B ^ Ktt 
decays given in Eq. (|l|) . It is clearly important to understand what we can learn 
from these results. Contributions from Fleischer, Gronau, and Neubert on this 
subject are rather technical. Nevertheless, it is an important technicality, as it 
must be understood when information is extracted from data. So, rather than 
summarizing what they have presented, I have presented necessary formalism 
which I hope is useful in following their work. 

Feynman graphs shown in Fig. |l| illustrates the class of operators which 
are generated by QCD and electroweak radiative corrections. The weak Hamil- 
tonian which causes these decays can be written as 

■H = ^|e„[Ci(M)or + C2{^im + ec[Ci(M)OJ + C2(m)0^] 

10 . 

-Ctl^aCA^)©, [+h.c., (10) 

where ^ ^*qb^qs, 0\ = (bs)v-A{m)v-A and Of = (bq)v-A{qs)v-A, 
Og = ^{bs)v-Aj2qeq{qq)v-A, OiQ = ^{bq)v~Aj2q^qi^^)v-A- Let us try 
to understand the isospin structure of these operators. The up tree graph Fig. 
|l|(a) contains u and u quarks and they generate both AI = 0, and AI = 1 
terms in the effective Hamiltonian. Fig. |l|(b), the charm tree graph, contains 
all isosinglet quarks and thus they generate AI = operator. Fig. |l|(c), the 
penguin, gives contribution which is proportional to X]g=M d s c Q^'^Q and it 
gives only AI = operator. Finally, Fig. |l|(d), the electroweak penguin, gives 
both AI = 0, and AI — 1 operators. 

Now we consider the isospin properties of the up tree, the operator which 
is generated by Fig. |l|(a). Because it contains both AI — 0, and AI — 1 
components, we want to separate the operator into two parts: 

+ c+{fi)[oi + o^] + C-{ti)[oi + oi] (11) 

where C± = C2 ± Ci and 0± ~ 5(^2 ± Oi). Then the first two terms on the 
right hand side cause AI = 1 transition and the last two terms cause AI = 
transition. 

Next we showJ-tiat the electroweak penguins, can be expressed interms of 
existing operatorsBt3. Note that the standard model predicts that Cj^si^b) 
are very small and they are negligible compared to C9,io(to6). The operators 
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with dominant coefRcients Og and Oio, for q ^ u, d, can be written as linear 
combinations of Oi'"^ and Oj ^' respectively: 



2 



(12) 



where Tij denotes the Hamiltonian which transforms as isospin 1. The operators 
above are defined as: 

o±^^{oi-oi), d± = l{oi + oi), 

Ol = l{0-2±01), = Cio ± Cg. 

In studying B — > Kit and tttt decays, it is important to classify final states 
in terms of strong interaction eigenstates, i.e. isospin states. This will allow 
us to take in to account of all rcscattering effects which have been discussed 
extensively in the literature. 

A{B+ ^7r+K")=B.+A.+A. =P + A-\PI,^ 
-V2A{B+ -> Tr°K+) =Bi+Ai-2A3=P + T + C + A+ lP%y^ + Pew 
-A{B^ ^ TT-K+) ^ bI - aI - As^ =P + T+ \P%y^ 
V2A{B° 7r"K°) = bI - aI+2A3 = P -C - |P|^ - Pew 

where Ai = {{K7r)i\ni\B) and Bi = {{Ktt) i\no\B) . We have also given the 
decay amplitudes in terms of amplitudes classified by Feynman graph struc- 
ture: tree graph (T), color suppressed tree graph (C), annihilation graph (A), 
penguin graph (P), electroweak penguin graph (Pew), and color suppressed 
electroweak penguin graph {Pew)- 

These decays together with their charge conjugate version constitue 8 phys- 
ical observables. Unlike K tttt decays Watson's theorem cannot be applied, 
and we cannot say rrmch about final state interaction phases for these ampli- 
tudes. We thus writeB: 



+ 5^e«(ft+03) i = 1,2,3. (13) 



Here we separate the contributions which are proportional to = l^ule*"^^ from 
those proportional to and ^t. at and Pi are final state strong interaction 



phases. It is trivial to write and bi in terms of matrix elements of the effective 
Hamiltonian Eq. (p2|). Note that there are 12 independent parameters in Eq. 
(|l3|), and only 8 decay rates for B Kn. 

In terms of matrix elements of the Hamiltonian, we have 

Bi=P; A^^C+Ql + C^Q^ Ai^C+Ql + C^Qt (14) 

where 

C± = ^uC± - = |eu|C±(e* - S±), (15) 

and Qj. is an appropriate matrix element {{Ktt)i\0±\B). We also record 
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P = C+Ql + c_g^ + uc+Ql + + ^tJ2 ^^Q^ (16) 

where 

C± = |^«|C±(e^^^-i(5±), Qi^{{K7T)i\d±\B), and ^ {{K7r)i\dl\B). 

So far, we have been quite general. Now, we shall go on to discuss the the 
hadronic matrix elements. What do we know about these matrix elements? 
Over the years, we have learned quite a bit. In particular, we have learned that 
classifying operators in terms of their topology. A, C, T, P, etc., gives us fairly 
accurate intuition in guessing the size of the matrix elements. For example, 
we guess that A, the annihilation graph, should be quite small compared to 
T because it is suppressed by the small probability that the spectator quark 
and b quark come within the range so that they can annihilate. Similarly, C 
is suppressed compared to T because of the color factor. These statements 
imply definite relationships between hadronic matrix elements which appear 
in Eq. (0). These relations should be checked experimentally. But, for 
the time being, we shall proceed and ask if these conventional wisdom would 
allow us to determine (j)3: the KM phase. The first simplification is A{B^ 
Tr+K°) = Bx + Ai + Ri \P\e^'l"' if A and Pew is negligible compared to 
P. The second simplification is that Q- transforms like a AI — operator 

_ 3 

in the limit of U spin symmetry. So, Qt vanishes in the SU(3) limit and is 
proportional to SU(3) breaking interaction. For our purpose, we neglect it. 
Then -3^3 = iPlea/ze^-^a/s (e^-^a _ 5^)^ where |P|e3/2 = -3|C„|C+q|. These 
considerations make analysis of _B+ K^t:^ simple: 

- V2A{B+ K+n°) = A{B+ K"tt+) - 3A, (17) 
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Neubert considers 



where A(/) — (pa — (pp. Fleischer considerslli 



(18) 



^ _ [Br(i?o ^ X+TT-) + Br(i3° ^ K-tt+)] ^^^^ 
Bt{B+ -> K^TT+) + Bt{B- ifn-) ' 

The decay K+ TT is bit more comphcated because we have to confront 

the contribution from Ai. It involves two complex amplitudes. He suplements 
the complexity by also considering 

Br(i?° ^ K+n-) - Bt(b" ^ K-n+)] 

^0 — =n 

Br(B+ ^ K"tt+) + Br{B- K\-) 

Detailed numerical analysis indicates that both of these methods are useful for 
determining <p^. In this discussion, I had to simplify the problem in order to 
present the essence. I refer the reader to the original contribution for complete 
analysis. It is clear that their contributions lead to much progress but much 
more work is necessary along this direction. For example, only subset of B — > 
Kn has been considered. There are 8 of them altogether! 

5 SUSY in B and K decays 

Predictions of the minimal supersymmetric theories (MSSM) is essentially 
same as those of the SM. If nature has chosen MSSM, we will not learn anything 
new from experiments on B and K decays. We should not be too discouraged 
by this though, as it is likely that she has chosen a theory which is more ele- 
gant than the MSSM. But, as long as we can not specify which theory nature 
has chosen, it is not an easy task to analyze its predictions. There are as 
many as 124 parameters in a non- minimal SUSY, and perhaps even more. 
Because B and K decays give stringent restrictions on flavor changing neutral 
current strengths, we shall focus general predictions of FCNC processes of a 
non-minimal SUSY - mostly penguin effects. 

In SUSY, there is a bosonic partner for each helicity of a quark. Here we 
begin by examining a 6 x 6 squark mass matrix of the MSSM. 



M 



~ U*™3/2(1 + ■^MlMv)M^^ ml'^^M^k^p 



(21) 
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where 




y^DLR = [\A\m„, + M*^ j MT- (22) 

1 is a 2 X 2 unit matrix, M.d,{u) is a mass matrix of a D{U) quark, and 
is a corresponding diagonalized matrix. Note that there new FCNC effects 
from additional flavor mixing among squarks. Others are standard MSSM 
parameters. To go from MSSM to a more general theory, let us identify 

AIj, = Amy^Mo (l + -^MIjMu 

Al^^A*my,(^l + ^MljMu^M^j,. (23) 

and generalize Sll = Slr = 6rl = ^ , and Srr = where 

TO is an average squark mass, as new arbitrary dimensionless parameters. We 
then study experimental constraints on 6 ignoring all theoretical prejudice. 

This analysis has been discussed in detail by Masiero. Bounds on 5 has 
been obtained from experimental information on various FCNC processes. For 
an average squark mass and gluino mass of 500GeV, the bounds on 5 ranges 
from 10^^ to 10^^. It sould be noted that the neutron edm gives a bound of 
Im {Slr)ii ~ 10~^ and it is natural to assume that other components of S are 
of the same order of magnitude. If this is the case, it may be difhcult of SUSY 
effects to show up in B and K decays. 

6 B decays to riK{K*), r]'K(K*), and r]'Xs 
CLEO has observedlli 

Br(B ^ 7/K+) ^ (^6.5^1'^ ± 0.9^ x 10"^ 
Br{B -> r]'K°) = (4-7^2 q ^ ^-^j ^ ^^^^ ' 



(24) 
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These branching ratios are surprisingly large. Du has presented a review of 
work in progress to undersand these large branching ratios. It is likely that 
these branching ratio is large because of gluonic content of 77'. 

Among various suggestions, a particularly interesting one is that of Soni 
and AtwQod. They attempt to compute 77' — > glue glue by considering triangle 
anomaly tJ. They then estimate b ^ rj' gluon. We should note, however, that 
major contribution to this decay comes from off shell gluon. Thus the validity 
of the anomaly calculation is questionable at best. A universal characteristic of 
all the work presented by Du is that each author picks up their favorite diagram 
and estimates its contribution. Nature does not work this way. They have to 
consider all possible diagrams and sum them up. Clearly global analysis is 
urgently needed. Also, there are large amount of data on the gluonic content 
of if. Such global analysis must be consistent with the previously known 
properties of 77'. 

7 A new calculation for direct CP violation: — 

Ever since the discovery of CP violation, experimentalists have been looking 
for an evidence of direct CP violation. Now that the result of NA31 has 
been confirmed by E832, the direct CP violation has been experimentally 
established. 

The challenge for theorists is to extract physics from the new value of ^. 
This is not an easy task. Before we compute CP violating amplitudes for 
K — > TTTT decay, we have to demonstrate that we understand CP conserving 
decay amplitudes. This means that we need to understand the A/ = i rule. 
At the moment there is no clear understanding_of this rule. So, one choice is 
to obtain hadronic matrix elements from datall3. If the AI — ^ rule is due 
to some new physics, this procedure will miss the new physics. Clearly, this is 
not satisfactory. We want to compute everything from basic principles. The 
approach taken by Paschos is an attempt along this direction. 

Let us start from the defining equation: 

,' = ^.(^.-^0) 1 ("^E.M. (25) 
2^/2 V^o A2J 

where Aj — {{TnT)j\H\K) and lu — \A2/Ao\ ^ .05 and Sj is the tttt phase shift 
for isospin / channel. 

In terms of operators in the effective Hamiltonian Ti, 

i = _,^KS.-^o) ^ 10-4 (^^) r y,{Q,), - 1 ^ y,(Q,)2 j (26) 
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where r = ^^^^^'^ = 336 GeV ^; j/j is the imaginary part of Wilson coeffi- 
cients; {Qi)j is a hadronic matrix element ((7r7r)/|0i|ii'); r = — . 

In this workshop, Paschos described computation of hadronic matrix ele- 
ments {Qi)i based on the -j^ expansion. They have obtained 

5 X 10"^ < - < 22 X 10"^ (27) 

for nis = (150 ~ 175)MeV. Their prediction increases if smaller values of 
is taken. 

This is an on going study and much more work remains: (I) The Wilson 
coefficients can be computed reliably only at some large energy scale Ac. But 
hadronic matrix elements can be computed only at low energy scale. So, 
there has to be some compromise. Thciv avc chosen Ac ~ SOOMeV. Wilson 
coefficient functions have large scale dependence in this region. When the 
coefficient functions and hadronic matrix elements are combined, the result for 
^ should not have the scale dependence. This has to be studied carefully. (II) 
It is necessary to demonstrate that K ^ tttt decay can be understood in this 
framework. Indeed, their result for — > tt+tt*^ decay, the AJ = | channel, 
is consistent with experiment. But it is necessary to understand the AI = ^ 
amplitude. Personally, I am skeptical toward a claim that the A7 = ^ rule can 
be understood within the frame work of expansion. 

8 Summary 

We have tried to have extended discussions on penguins. We tried to under- 
stand how we might extract </)2 and 03 from data. We don't think there is one 
best method to extract these angles. It is an experiment driven field, and time 
will tell. We tried to understand — from basic principles of the SM. There 
is much more work to be done along this direction. We tried to understand 
B ^ T]' + X decays. This also requires more work. Seeing effects of SUSY in 
B and K decays is an exciting possibility. We have to keep on searching. 
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